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Preface 


The  Liquid  Droplet  Radiator  CLDR)  consists  of  a  column 
or  sheet  of  liquid  droplets  moving  through  space  from  a 
droplet  generator  to  a  collector.  The  droplets  carry  the 
waste  heat  generated  by  a  space  power  system  and  radiate 
this  waste  heat  directly  to  space  during  their  flight.  The 
liquid  droplets  are  collected  at  a  lower  temperature, 
reheated,  and  pumped  to  the  generator  and  reused  to  continue 
to  remove  waste  heat  from  the  thermodynamic  power  cycle. 

This  study  was  a  parametric  analysis  of  a  cylindrical  LDR  to 
estimate  its  performance  and  operating  characteristics  using 
a  varying  pump  specific  mass. 

I  offer  a  well  deserved  expression  of  gratitude  to  my 
advisor,  Lt  Col  Ronald  Tuttle  of  the  Air  Force  Institute  of 
Technology  at  Wright-Patterson  AFB,  Ohio,  for  his 
instrumental  guidance  and  patience  during  this  project. 

Also,  I  would  like  to  thank  Dr.  Alden  Presler,  Dr.  Robert 
Siegel,  Ms.  Carolyn  Coles,  and  Mr.  Alan  White  of  the  NASA 
Lewis  Research  Center,  Cleveland,  Ohio,  for  their  assistance 
and  for  providing  current  research  information  on  the  LDR. 
Ultimately,  I  must  give  my  most  heartfelt  appreciation  to  my 
wife,  Linda,  for  her  understanding  and  her  ability  to  manage 
our  family  and  three  children  with  little  help  from  me  for 
the  past  year  and  a  half. 


Gerald  L.  Buckner 
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Abstract 


JThis  study  was  a  parametric  investigation  of  the 


performance  and  operating  characteristics  of  a  cylindrical 


Liquid  Droplet  Radiator  CLDR)  for  use  in  space.  The  LDR 


system  mass  per  heat  radiated  was  minimized  as  a  function  of 


the  average  droplet  temperature  at  the  collector.  This 


study  was  similar  to  the  work  of  Karl  Knapp  C1980),  however 


a  new  pump  specific  mass  term  wa3  used  in  the  total  system 


mass  calculation.  Knapp  used  a  pump  specific  mass  defined 
as  kg-sec*' 3 /nr  .  'This  study  used  a  pump  specific  mass  term 


defined  as  pump  mass  per  liquid  mass  flow  rate  to  develop  a 


physically  meaningful  pump  specific  mass  term  for  use  by 


design  engineers.  '  The  new  pump  specific  mass  was  varied 


from  10.0  kg/ (kg/ sec)  to  40.0  kg/ C kg/ sec),  based  on 


available  industry  standard  pumps.,-  The  average  droplet 


temperature  at  the  collector  was  calculated  to  minimize  the 

KiLooutft  KiLov/trj 

LDR  system  mass  for  heat  loss  rates  of  10  Kw,  100  Kw,  1000 
Kw,  and  10,000  Kw  for  mission  lifetimes  of  zero,  ten, 


twenty,  and  thirty  years.  The  initial  droplet  temperature 


at  the  generator  was  fixed  at  300  degrees  Kelvin.  A  silicon 


oil,  Trimethylpentaphenyltrisiloxane  CDOW  705),  was  modeled 


due  to  its  low  vapor  pressure  C approx.  3  x  10  8  Pa)  at  300 


degrees  Kelvin.  The  variable  pump  specific  mass  term  offers 


the  design  engineer  a  range  of  possible  pump  masses 


depending  on  system  pressure  and  flow  rate  requirements. 
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THE  LIQUID  DROPLET  RADIATOR  IN  SPACE:  A  PARAMETRIC  APPROACH 

I.  INTRODUCTION 

A  successful  space  mission  must  have  a  source  of 
electrical  power  whether  the  mission  is  manned,  unmanned, 
scientific,  or  nationally  strategic.  The  generation  of  this 
electric  power  will  require  the  rejection  of  waste  heat. 

For  example,  the  Strategic  Defense  Initiative  will  h&ve 
space  based  systems  generating  large  amounts  of  electrical 
energy  with  much  waste  heat  energy  to  be  radiated  to  space. 
Other  space  applications  requiring  from  10.0  Kilowatts  to 
over  100  megawatts  include:  Space  Based  Radars, 
Nuclear/Electric  Orbital  Transfer  Vehicles,  Space  Based 
Weapon  Systems,  and  the  Space  Station To -operate 
these  space  power  systems,  a  lightweight  highly . survivable 
radiator  is  considered  the  "enabling  technology"  (2:9.1).  A 
Liquid  Droplet  Radiator  (LDR)  can  be  seven  times  lighter 
than  conventional  heat  pipe  radiators  of  similar  size 
(4:14).  Thus  the  LDR  technology  development- can  enable  the 
Air  Force  and  NASA  to  conduct  their  space  missions  with  less 
mass  in  orbit  and  with  resulting  cost  savings. 


Background 


A  space  power  system  must  operate  with  the  rejection  of 


waste  heat  from  the  thermodynamic  cycle  used  to  generate  the 
electrical  energy  for  use  in  space.  To  reject. this  waste 
heat  in  space,  radiators  must  be  designed  to  operate 


i 


I 

I 


efficiently  in  the  .harsh  space  environment.  Traditional 
radiator  designs  for  operation  in  space  are  the  pure  thermal 
radiator  and  the  bimodal  radiator.  The  thermal  radiator 
Cmetallic  fins  for  example)  has  been  designed  in  various 
sizes  to  maximize  area,  minimize  weight  and  size,  and 
maximize  protection  from  meteoroid  damage.  In  the  bimodal 
radiator,  a  small  thermal  radiator  is  used  to  reject  heat 
for  low  power  operations  while  for  high  power  operations,  a 
dump  fluid,  such  as  hydrogen,  is  used  to  remove  excess  heat 
by  dumping  it  directly  into  space.  These  traditional 
radiators  are  not  acceptable  to  meet  the  heat  rejection 
requirements  of  multimegawatt  space  power  systems  of  the 
future  since  radiator  mass  requirements  will  dominate  these 
large  power  systems  making  them  prohibitively  massive  and 
very  expensive.  Less  massive  radiator  designs  are  needed  to 
remove  the  increased  waste  heat  from  multimegawatt  power 
systems  for  efficient  space  power  generation.  The  liquid 
droplet  radiator  is  one  possible  proposed  design  C3:i).  The 
radiator  coolant  liquid  is  formed  into  aerosol  droplets  and 
passed  as  a  column  or  thin  sheet  of  droplets  from  a  droplet 
generator  through  the  vacuum  of  space  to  a  collector  CFigure 
1).  The  waste  heat  is  rejected  to  space  by  radiation  from 
each  droplet  while  in  flight.  The  advantages  of  such  a 
design  include  large  droplet  surface  areas  per  mass, 
lightweight  fluid  mass,  and  lower  sensitivity  to  meteoroid 
damage  (4:1). 
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Figure  i.  S  thematic  of  Two  Liquid  Droplet  Radiators 
(adapted  from  References  1  and  9) 


Obj  ective 


T 


he  ooj ective  of  this  study  was  to  investigate  the 


performance  and  operating  characteristics  of  a’ cylindrical 
LDR  for  use  in  space  by  minimizing  the  mass  per  heat 
radiated  as  e  function  of  the  average  droplet  temperature  at 
the  collector  using  a  new  pump  specific  mass  terifi  defined  as 
pump  mass  per  liquid  mass  flow  rate.  ^The^new  pump  specific 
mass  was  varied  from  10.0  kg/Ckg/sec)  to  40.0  kg/Ckg/sec), 
bai.ed  on  available  industry  standard  pumps,  to  offer  the 


design  engineer  a  range  of  possible  pump  masses  depending  on 
system  pressure  and  flow  rate  requirements. 


A  silicon  oil,  Trimethy lpentaphenyl Lrisiloxane  (DOW 

705),  was  used  in  this  study  in  a  low  temperature  heat 

rejection  application  with  a  300  degrees  Kelvin  droplet 
»  , 

temperature  at  the  generator.  The  LDR  stream  length,  stream 
diameter,  droplet  radius,  and  total  system  mass  per  radiated 
power  were  found  for  the  minimum  mass  LDR  system  as  a 
function  of  the  average  droplet  temperature  at  the  collector 
for  heat  rejection  rates  of  10  Kw,  100  Kw,  1000  Kw,  and 
10000  Kw  with  mission  operation  times  of  zero,  ten,  twenty, 
and  thirty  years. 

Approach 

The  total  LDR  system  specific  mass  was  calculated  as  a 
function  of  the  average  temperature  at  the  collector  for  a 
given  set  of  fixed  parameters  and  for  varying  pump  specific 
mass  values.  The  temperature  corresponding  to  the  minimum 
mass  system  was  then  used  to  calculate  droplet  stream 
length,  diameter,  and  droplet  radius  for  each  case. 

Overview 

The  theory  of  operation  of  the  LDR  droplet  column  is 
presented  in  Chapter  II.  The  parameter  investigation  is 
presented  in  Chapter  III  for  the  new  pump  specific  mass 
term.  The  results  are  presented  in  Chapter  IV  and 
conclusions  and  recommendations  are  presented  in  Chapter  V. 


IX.  THEORY 


The  LDR  in  this  study  was  modeled  as  a  cylindrical 
stream  of  droplets  (Figure  1).  The  procedural  approach  for 
this  study  was  adapted  from  Reference  1.  The  radiating 
stream  is  assumed  to  be  shadowed  so  solar  absorption  need 
not  be  considered  (1:3).  A  single  droplet  radiates  heat  as 
it  travels  through  space  and  at  any  time  this  heat  loss  is 
given  by: 


q 


(47ra2)aFT4 


(1) 


where 


4  =  droplet  heat  loss  rate  to  space  (joules/second) 
a  =  droplet  radius  (meters)  2  i| 

a  =  Stefan-Boltzmann  Constant  (5.67  x  10“8  watts/m  ~K  ) 
F  =  average  gray  body  view  factor  for  droplet  at  stream 
center  (less  than  one) 

T  =  absolute  droplet  temperature  at  any  time  (Kelvin) 


This  equation  models  the  droplet  as  a  gray  body  with 
constant  average  emissivity.  The  instantaneous  radiation 
rate  is  equal  to  the  rate  of  energy  loss  resulting  in  this 
equation: 


2  4 

47Ta  <rFT 


dT 

dt 


(2) 


where 


C  ~  specific  heat  capacity  (1670  .1oules/kg  -  K) 
p  =  density  of  droplet  (1000  kg/m3 ) 
t  =  droplet  transit  time  (seconds) 


The  integration  of  Eq  C2)  and  solving  for  time  results  in 
this  equation  for  the  droplet  transit  time: 


y*47ra2<7F  dt  =  si  ~ 

.  cpa  1 _ 1_ 

9*F  t3  t3 


where 


C33 


=  initial  droplet  temperature  CKelvin) 


The  transit  time  from  generator  to  collector  is  equal  to 
L/V  where  L  is  the  length  of  the  droplet  stream  in  meters 
and  V  is  the  droplet  velocity  in  meters/ second.  The 
temperature,  T,  is  equal  to  T2  at  the  collector  and  Eq  C3) 
can  be  solved  for  the  droplet  radius  Ca)  by  substituting 
t=L/V  and  T=T2  resulting  in  this  equation: 


a 


9<7FL 

CpV 


dyi^r3  -  1 


C4) 


The  average  rate  of  heat  loss  from  each  droplet  is  given  by: 


/•L/V 

/ 

X  J  o _ 

^  T. 


L 

V 


2  4  3(1-1,/!.) 

=  4ma  aFT,  - -  x 

-L  /  m  /m  \ 


C5) 


(t2/Ti) 


-  1 


For  the  stream  volume  with  the  number  of  droplets  per  cubic 
meter  equal  to  n  and  the  stream  diameter  given  by  D  in 
meters,  the  total  rate  of  heat  loss  for  the  stream  is  given 
by  the  following  equation: 
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*  2  2  2  4 

Q  =  n7Ta  D  L<rFT* 


1 


ad-Tj/T^ 


<  W 


The  gray  body  view  factor  from  Reference  1  is  given  as: 


i+  ±-  1 

V  € 


where 


e  a  the  em.  ;3ivity  of  the  droplet  surface 

*  the  conservative  black  body  view  factor  for  the 
droplet  at  stream  center 


The  value  for  eta  is  used  as  a  conservative  average  value 


for  all  the  droplets  and  is  given  as  Cl: 4): 


*7  =  1- 


2  ,  2 
n  nDa 


Solving  Eq  C8)  for  n  and  substituting  with  Eq  C7)  into  Eq 


C6),  the  expression  for  the  total  heat  loss  from  the  droplet 


stream  is  found  to  be: 


4  3(1-T  /T  > 

Q  *  BDhaT?  - — - - 

(Ty'T  )“3  -  1 


where 


1  .  1 


~  +  ~  -  1 
n  e 


The  value  for  B  is  substituted  into  Eq  C9)  and  the 


tv 


resulting  equation  is  differentiated  with  respect  to  eta  and 
the  result  set  equal  to  zero  in  order  to  reasonably  optimize 
the  heat  flow  per  unit  area  resulting  in  this  equation  for 
eta: 


The  total  mass  of  the  system  results  from  the  mass  of  the 
liquid  to  carry  the  waste  heat,  the  mass  of  the  generator 
and  collector,  the  pump  mass,  and  the  extra  liquid  mass 
needed  to  replace  evaporation  losses.  The  total  mass  of  the 
system  is  given  by  (adapted  from  Reference  1): 


"r 


..  Q(1+B)L  .  ir(mc)Q2 

a  cv^f^T"!  +  -"IS — 


(T2/T1)"3  -  1 
BLaT3(T1-T2) 


+ 


Q(mp) 

c(t1-t2) 


+ 


nQ  P(T1)(T1/Tv)t 
FaT3  (T1-T2)(2ttRT1)1/2 


(11) 


where 


3  =  parameter  variable,  ratio  of  liquid  mass  in 
reheating  station  to  mess  in  stream 
me  =  parameter  variable,  specific  mass  of  drople£ 
generator  and  collector  per  unit  area  (kg/m  ) 
mp  =  parameter  variable,  pump  specific  mass 
(kg/ (kg/ sec) ) 

P(T^)  =  vapor  pressure  of  droplets  at  T.,  =  300  Kelvin 
(pascals)  1 

Tv  =  value  used  to  fit  vapor  pressure  curve 
C 15000  Kelvin) 

R  =  gas  constant  (15.22  j oules/kg-Kelvin) 
x  =  mission  time  (years) 


The  first  term  gives  the  mass  of  liquid  needed  to  carry  the 
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amount  of  waste  heat  in  the  stream  and  the. fraction  of  the 
total  liquid  mass  in  the  reheating  station.  The  second  term 
is  the  mass  of  the  droplet  generator  and  collector.  The 
third  term  is  for  the  pump  mass  and  includes  the  new  pump 
specific  mass  Cmp)  in  terms  of  pump  mass  per  mass  flow  rate. 
The  fourth  term-  is  the  total  mass  lost  due  to  evaporation  or 
the  amount  of  additional  liquid  mass  needed  to  replace  that 
lost  to  evaporation  Cl: 8). 

To  find  the  stream  length  L  to  minimize  the  total 
mass  for  a  given  set  of  the  other  parameters,  Eq  C 11)  is 
differentiated  with  respect  to  L  and  set  equal  to  zero 
resulting  in  this  equation: 


ttCVTjQ  [(T2/T;i)~3  -  lj2(mc) 
18B2( 1+B ) (otJ) 2( 1-T2/T1) 


(12) 


This  equation  for  L  is  substituted  into  Eq  C 11)  and  when 
the  total  heat  rejection  rate  is  divided  through  the 
equation,  the  following  equation  for  the  specific  mass  of 
the  system  in  terms  of  kg/Kw  is  found  (adapted  from 
Reference  1): 


JJt  _  3Q1/3 

Q  2(1-T2/T1) 11/3 

(mp)  nPCT^CyTvlx 

C(T1:V  FoT  J(  1-T2/T1)  (  2tiRT1)  1/2  (13) 
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Possible  inaccuracies  in  the  direction  of  flight  of  the 
liquid  droplets  may  cause  the  stream  to  widen  at  the 
collector  and  droplet  mass  will  be  lost.  To  minimize  this 
stream  widening,  the  length  to  diameter  ratio  (L/D)  is  held 
fixed  at  a  recommended  value  of  250  (1:10).  Equation  Cll) 
is  rewritten  by  using  Eq  (9)  resulting  in  Eq  (14)  below 
where  the  only  dimensions  appearing  are  in  the  ratio  L/D 
(adapted  from  Reference  1): 


(1+3) (L/D) 1/2 

(T2/T1)“3  -  l]1/2  Q1/2 

CVT1( 3BaT^ 

)1/2(i-t2/t1)3/2 

+  (mp) 

+  6(t“-t2) 


•n(mc) 

(T2/Ti)-3  .  i] 

12BotJ 

(l/d)(i-t2/t1) 

nPCTj.MTj/Tv)! 

FotJ(2ttRT1)1/2(1-T2/T1) 


(14) 


The  stream  diameter  is  found  from  Eq  (9)  for  a  fixed  L/D 
ratio  as: 


Q1/2[(t2/t1)“3  -  i]1/2 
|^3B(l/d)otJJ1/2(i-t2/t1)1/2  C15‘ 

The  length  L  of  the  droplet  stream  is  found  for  the  same 

fixed  L/D  ratio  by  multiplying  the  value  of  D  from  Eq 

(15)  *imes  L/D.  The  droplet  radius  a  is  found  from  Eq 

(4)  after  L  is  calculated. 

The  stream  velocity  is  found  from  an  equation  for  the 
ratio  of  droplet  kinetic  energy  to  the  total  heat  rejection 

$ 
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rate  Cl: 6)  as: 


Y  2C(T1-T2)  C16) 

or  solving  for  the  stream  velocity 

V  .  (2CY(T1-T2))1''2  C17) 

The  value  of  gamma  is  required  to  be  much  less  than  one  to 
reduce  this  ratio  and  thus  keep  the  pump  power  within  a 
small  fraction  of  the  total  heat  rejection  rate  Cl: 6). 

The  equations  in  this  chapter  represent  LDR  operation 
under  steady-state  conditions.  A  detailed  design  of  a  LDR 
would  involve  many  other  factors  such  as  the  wear  of  moving 
parts,  problems  of  corrosion  and  erosion,  deployment  and 
start  up  techniques,  the  effects  of  variable  heat  loads, 
temporary  shutdowns,  and  possible  on-orbit  servicing  Cl: 10). 


III.  PARAMETER  INVESTIGATION 

Parameter  values  were  selected  for  use  in  Eq  CIS)  and 
the  total  system  specific  mass  was  calculated  as  a  function 
of  the  average  droplet  temperature,  T2,  at  the  collector  for 
various  pump  specific  mass  values.  The  ratio,  L/D,  was 
calculated  using  the  droplet  temperatures  corresponding  to 
the  minimum  mass  system.  If  the  L/D  value  was  greater 
than  250,  a  new  minimum  specific  mass  was  calculated  using 
Eq  C14)  with  L/D  fixed.  The  parameter  values  used  for  the 
calculations  are  shown  in  Table  I  and  were  adapted  from 
Reference  1  except  for  pump  specific  mass  Cmp).  Refer  to 
Appendix  B  for  calculation  of  some  typical  mp  values. 

TABLE  I 

Parameters  Used  for  Calculations 


*  optimistic  value  **  realistic  value  (1:13) 
12 
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The  total  system  specific  mass  was  found  using  Eq  C13) 
or  Eq  C14)  for  fixed  L/D  =  250  ratio  for  a  range  of  pump 
specific  mass  values  as  shown  in  Table  I  for  the  two  values 
of  3  and  me  using  the  computer  program  in  Appendix  A. 

The  curves  generated  by  these  calculations  are  shown  in 
Figures  2-17  for  heat  loss  rates  of  10  Kw,  100  Kw,  1000  Kw, 
and  10000  Kw  and  mission  times  of  zero,  ten,  twenty,  and 
thirty  years.  The  parameter  values  used  in  Figures  2-17 

p 

were  3  =0.10  and  me  =  40.0  kg/m  equivalent  to  5 

millimeters  of  steel  and  called  optimistic  values  C1:A-1). 

A  similar  set  of  data  Cno  figures)  was  calculated  with  3  = 

2 

0.20  and  me  =  100.0  kg/m  equivalent  to  12  millimeters  of 
steel  and  called  realistic  values  Cl:A-3).  In  both  cases 
y  =  0.005  and  e  =0.7  for  the  calculations.  No  figures 
were  generated  for  the  second  set  of  data  to  avoid 
redundancy.  The  temperature  values  tabulated  in  Chapter  IV 
for  T2  were  those  values  corresponding  to  the  minimum 
points  on  each  of  the  curves.  These  temperatures  are  the 
average  droplet  temperatures  at  the  collector  after  the 
droplets  radiated  their  waste  heat  to  space  as  a  function  of 
the  pump  specific  mass  for  the  minimum  mass  system.  The 
minimum  point  on  each  of  the  curves  in  Figures  2-17  shows  a 
reduction  in  the  droplet  temperature  at  the  collector  with 
increasing  pump  specific  mass  and  increasing  total  system 
specific  mass. 


A  comparison  of  the  average  droplet  temperature,  Tg,  at 
the  collector  for  the  minimum  mass  case  as  a  function  of 
pump  specific  mass  is  shewn  in  Figures  18-21  for  each  3 
and  me  parameter  value  for  various  mission  times  and  heat 
loss  rates.  In  Figures  18-21  the  lower  four  curves  were 
obtained  using  the  optimistic  values  of  3  and  me  from 
Table  I.  The  upper  four  curves  were  obtained  using  the 
realistic  values  of  3  and  me  .  The  symbols  on  the  curves 
are  the  actual  temperature  data  points  found  as  a  function 
of  pump  specific  mass  for  the  minimum  mass  system  for  the 
indicated  times  and  heat  loss  rates. 

In  Figures  22-25  a  comparison  of  the  system  specific 
mass  Ckg/Kw)  as  a  function  of  pump  specific  mass,  mp, 
Ckg/Ckg/sec) )  is  shown  for  the  indicated  heat  loss  rates  and 
mission  times.  In  Figures  22-25  the  lower  four  curves  were 
obtained  using  the  optimistic  values  of  fl  and  me  from 
Table  I,  while  the  upper  four  curves  were  found  using  the 
realistic  values  of  3  and  me.  The  symbols  on  the  curves 
are  the  actual  specific  mass  data  points  found  as  a  function 
of  pump  specific  mass  for  the  minimum  mass  system. 

The  system  specific  mass  as  a  function  of  the  heat  loss 
rate  is  shown  in  Figure  26  for  the  indicated  mission  times. 
Again,  the  lower  four  curves  were  obtained  using  the 
optimistic  values  of  3  and  me  from  Table  I,  while  the 
upper  four  curves  were  found  using  the  realistic  values  of 
3  and  me  .  This  figure  is  for  a  pump  specific  mass  value 
of  1C.0  kg/Ckg/sec). 
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The  following  key  is  for  the  symbols  on  the  curves  in 
Figures  2-17  representing  the  different  pump  specific  mass 
values  used  to  generate  the  curves. 


Symbol 

mp  Value 

□ 

10.0 

A 

15.0 

o 

20.0 

* 

25.0 

* 

30.0 

+ 

35.0 

X 

40.0 
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Figure  3.  System  Mass  per  Heat  Loss  vs  Droplet  Temperature: 
CIO  Kw,  Time  =  10  years) 


HEAT  LOSS 


TIME  - .20  YEARS 


200.00 

220.00 

240.00  260.00 

280.00 

DROPLET  TEMP2  (K) 

Figure  4.  System  Mass  per  Heat  Loss  vs  Droplet  Temperature 
CIO  Kw,  Time  =  20  years) 
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Figure  5. 


System  Mass  per  Heat  Loss  vs  Droplet  Temperature: 
CIO  Kw,  Time  =  30  years) 


r 
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Figure  7.  System  Mass  per  Heat  Loss  vs  Droplet  Temperature 
(100  Kw,  Time  =  10  years) 


Figure  8. 


System  Mass  per  Heat  Loss  vs  Droplet  Temperature: 
C100  Kw,  Time  =  20  years) 
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Figure  10. 
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System  Mass  per  Heat  Loss  vs  Droplet  Temperature: 
(1000  Kw,  Time  =  0) 
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Figure  19. 


Droplet  Temperature  vs  Pump  Specific  Mass: 
(100  Kw) 


PI 


Figure  20.  Droplet  Temperature  vs  Pump  Specific  Mass 
< 1000  Kw) 
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IV.  RESULTS 


The  numerical  results  are  shown  in  Tables  II-IX  with 
Knapp's  values  also  listed  Cl).  The  average  droplet 
temperature  at  the  collector,  total  LDR  specific  mass,  LDR 
length,  diameter,  and  droplet  radius  were  calculated  for 
each  value  of  the  new  pump  specific  mass  Cmp)  for  the 
minimum  mass  system  for  four  different  heat  loss  rates  and 
four  different  mission  times. 

10  Kilowatt  Heat  Loss 

In  the  optimistic  case,  Table  II,  Knapp's  values  were 
within  four  percent  of  the  LDR  temperature,  length,  and 
diameter  values  for  the  10.0  and  20.0  pump  specific  mass 
Cmp)  values.  However,  Knapp '3  results  for  total  system 
specific  mass  best  fit  the  mid  range  of  mp  values  C20.0  to 
25.0).  For  the  droplet  radius,  Knapp's  values  were 
approximately  one  half  the  radius  values  for  the  low  mp 
va.w.es. 

The  opposite  occurred  using  the  realistic  parameter 
values,  Table  III,  where  Knapp's  results  for  temperature, 
length,  and  diameter  values  agreed  within  two  percent  for 
the  mp  value  of  40.0,  while  the  total  system  specific  mass 
was  within  thirteen  percent  agreement.  Knapp's  droplet 
radius  values  were  again  approximately  one  half  the  radius 
values  for  an  mp  value  of  40.0.  There  was  an  apparent 
typing  error  in  Knapp's  report  for  the  droplet  radius  for 
the  ten  year  mission  time. 
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The  results  in  Figure  18  show  a  lower  droplet 
temperature  at  the  collector  with  increasing  pump  specific 
mass,  while  Figure  22  results  indicate  a  corresponding 
increase  in  the  total  system  specific  mass  with  increasing 
pump  specific  mass. 

100  Kilowatt  Heat  Loss 

In  the  optimistic  case,  Table  IV,  Knapp's  values  for 
total  specific  mass,  temperature,  length,  and  diameter  best 
agreed  within  three  percent  for  the  lower  mp  value  of 
10.0,  except  for  the  droplet  radius.  Knapp's  droplet  radius 
was  approximately  one  half  the  droplet  radius  for  the  mp 
value  of  10.0. 

A  similar  result  was  observed  for  the  realistic  values, 
Table  V,  where  Knapp's  results  for  total  specific  mass, 
temperature,  length,  and  diameter  were  in  very  good 
agreement  within  five  percent  for  the  mp  values  of  15.0 
and  20.0.  For  this  range  of  mp  values,  Knapp's  droplet 
radius  was  low  again  by  approximately  one  half  the 
calculated  radius  values.  There  were  apparently  typing 
errors  in  the  Knapp  paper  for  the  stream  diameter  at  twenty 
and  thirty  year  mission  times. 

The  results  in  Figure  19  show  a  lower  droplet 
temperature  at  the  collector  with  increasing  pump  specific 
mass,  while  Figure  23  results  indicate  a  corresponding 
increase  in  the  total  system  specific  mass  with  increasing 
pump  specific  mass. 


1000  Kilowatt  Heat  Loss 


In  the  optimistic  case,  Table  VI,  Knapp's  values  for 
temperature,  total  specific  mass,  length,  and  stream 
diameter  were  found  to  agree  within  eight  percent  with  all 
values  for  a  mp  value  of  10.0,  except  for  the  droplet 
radius.  The  Knapp  radius  values  were  again  approximately 
one  half  the  calculated  radius  values  for  the  mp  value  of 
10.0. 

For  the  realistic  case,  Table  VII,  Knapp's  values  of 
temperature,  total  specific  mass,  length,  and  stream 
diameter  were  found  to  agree  within  two  percent  with  all 
values  for  a  mp  value  of  10.0,  except  for  the  droplet 
radius.  The  radius  was  again  approximately  one  half  the 
calculated  radius  values  for  the  mp  value  of  10.0 

The  results  in  Figure  20  show  a  lower  droplet 
temperature  at  the  collector  with  increasing  pump  specific 
mass,  while  Figure  24  results  indicate  a  corresponding 
increase  in  the  total  system  specific  mass  with  increasing 
pump  specific  mass. 

10000  Kilowatt  Heat  Loss 

For  the  optimistic  case,  Table  VIII,  Knapp's  values  for 
temperature,  total  specific  mass,  length,  and  stream 
diameter  were  found  to  agree  within  seven  percent  with  all 
values  for  a  mp  value  of  10.0,  except  for  the  droplet 
radius.  The  Knapp  study  droplet  values  were  approximately 
three  fifths  of  the  calculated  radius  values  for  the  mp 
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value  of  10.0. 


For  the  realistic  case,  Table  IX,  Knapp's  values  of 
temperature,  total  specific  mass,  length,  and  stream 
diameter  were  found  to  agree  within  three  percent  with  all 
values  for  a  mp  value  of  10.0,  except  for  the  droplet 
radius.  For  the  mp  value  of  10.0,  Knapp's  radius  values 
were  approximately  one  half  the  calculated  radius  values. 

The  results  in  Figure  21  show  a  lower  droplet 
temperature  at  the  collector  with  increasing  pump  specific 
mass,  while  Figure  25  results  indicate  a  corresponding 
increase  in  the  total  system  specific  mass  with  increasing 
pump  specific  mas3. 
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TABLE  II 


Results  for  10  Kw  Using  Optimistic  Values 


Mission 

Time 


.  Pump  Specific  Mass  C kg/ C kg/ sec)) _ 

10. Oj  15. 0 1  20.0  I  25. ol  30.01  35.01  40.0  Knapp 


T2=|  241 


10  years 


225  221  217  214 


a=  7.0  0. 5 


239 


20  years 


30  years 


Mt/Q=  0.387  0.435  0.479  0.520  0.559  0.596  0.631  0.501 
L=  71.8  73.8  75.9  77.4  78.9  80.6  81.8  70.8 


D=  0. 287  0.  295  0.  304  0. 310  0.  316  0. 322  0. 327 


4.1  3.9 


224  220  217  214 


0. 401  0.  448  0. 491  0. 531  0. 570  0.  606  0. 642  I  0. 511 


72.5 


77.8  79.3  80.6  81.8 


0.  290  0.  298  0. 305  0. 311  0. 317  0. 322  0. 327  I  0. 291 


6. 2  5.. 


232  227 


4.1  3.9 


6.0  5.3 


236  230 I  225  22l|  218|  215l  212 


213 

231 

0.652 

0.523 

82.2 

74.  8 

0.329 

0.299 

3.  8 

2.  8 

212 

228 

0.  662 

0.536 

82.7 

75.  9 

0.  331 

0.  303 

3.7 

2.7 

6.7  5.7  5.0  4.6  4.2  4.0  3.7 


Units:  T2  CK),  M<p/Q  Ckg/Kw),  L  Cm),  D  Cm),  a  Cmicrons) 
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TA  'LE  III 


Results  for  10  Kw  Using  Realistic  Values 
CBeta  =0.20,  me  =  100.0) 


Mission 

Pump  Specific  Mass 

C kg/ C kg/ sec) 

) 

Time 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

40.  0 

Knapp 

n 

CM 

H 

o 

257 

252 

248 

244 

241 

238 

236 

Mj/Q= 

0.  747 

0.  813 

0.  873 

0.  929 

0.  981 

1.031 

1.078 

L= 

67.  1 

68.  5 

69.  7 

70.  9 

71.  8 

72.  8 

73.  5 

D= 

0.  268 

0.  274 

0.  279 

0.  284 

0.  287 

0.  291 

0.  294 

a= 

13.  3 

11.  1 

9.  7 

8.  5 

7.  7 

7.  1 

6.  7 

3.  2 

10  years 

256 

251 

247 

243 

240 

237 

235 

234 

0.  767 

0.  831 

0.  890 

0.  944 

0.  996 

1.044 

1.092 

1.  23 

67.4 

68.  8 

70.0 

71.  2 

72.  2 

73.  2 

73.  8 

73.  8 

0.  270 

0.  275 

0.  280 

0.  285 

0.  289 

0.  293 

0.  295 

0.  295 

12.  8 

10.7 

9.  3 

8.2 

7.  5 

6.  9 

6.  5 

*8.0 

20  years 

254 

— 

250 

246 

242 

239 

237 

234 

234 

0.  787 

0.  849 

0.  906 

0.  960 

1.010 

1.059 

1.  105 

1.  24 

67.  9 

69.  1 

70.  3 

71.  5 

72.  5 

73.  2 

74.  2 

73.  8 

0.  272 

0.  276 

0.  281 

0.  286 

0.  290 

0.  293 

0.  297 

0. 295 

11.  9 

10.  3 

9.  0 

8.0 

7.  3 

6.  9 

6.  3 

3.0 

30  years 

253 

248 

245 

241 

239 

236 

233 

234 

0.  805 

0.  866 

0.  922 

0.  975 

1.025 

1.072 

1.  118 

1.  25 

68.  2 

69.  7 

70.  6 

71.  8 

72.  5 

73.  5 

74.  5 

73.  8 

0.  273 

0.  279 

0.  282 

0.  287 

0.  290 

0.  294 

0.  298 

0.  295 

11.  5 

9.  7 

8.  8 

7.  7 

7.3 

6.  7 

6.  2 

3.0 

Units:  Tg  CKO,  Mp/Q  (kg/Kw),  L  Cm),  D  Cm),  a  Cmicrons) 
*  possible  typo  in  Knapp  study 
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TABLE  IV 


t 


Results  for  100  Kw  Using  Optimistic  Values 
•CBeta  =  0.10,  me  =  40.0) 


Mission 

Pump  Specific 

Mass 

Ckg/Ckg/sec)) 

Time 

KHail 

mill 

*| 

my; 

mm 

Baai 

■Oil 

Knapp 

ii 

CM 

E- 

o 

233 

228 

224 

220 

217 

r 

214 

211 

237 

mt/q= 

0.  438 

0.481 

0.  522 

0.  560 

0.  597 

0.  632 

0.  667 

0.  437 

L= 

236 

241 

246 

251 

255 

259 

263 

230 

D= 

0.  943 

0.  965 

0.984 

1.004 

1.020 

1.030 

1.050 

0.  920 

a= 

19.  5 

17.  1 

15.  5 

14.  1 

13.  1 

12.  2 

11.  4 

10.0 

10  years 

232 

227 

223 

219 

216 

213 

211 

234 

0.  451 

0.  493 

0.  533 

0.  571 

0.  607 

0.642 

0.  676 

0.  450 

237 

242 

247 

252 

256 

260 

263 

233 

0.  947 

0.  970 

0.  989 

1.009 

1.024 

1.040 

1.051 

0.  933 

19.0 

16.  7 

15.  1 

13.7 

12.  8 

11.  9 

11.  4 

9.  5 

20  years 

230 

226 

222 

218 

215 

213 

210 

234 

0.  463 

0.  505 

0.  544 

0.  582 

0.  618 

0.  652 

0.  686 

0.  463 

239 

244 

248 

253 

257 

260 

264 

233 

0.956 

0.  974 

0.  994 

1.014 

1.030 

1.040 

1.057 

0.  933 

18.0 

16.  3 

14.  7 

13.  4 

12.  5 

12.0 

11.  2 

9.  5 

30  years 

229 

225 

221 

217 

214 

212 

209 

231 

0.  476 

0.  517 

0.  555 

0.  592 

0.628 

0.  662 

0.  696 

0.  476 

240 

245 

250 

255 

259 

261 

266 

237 

0.  960 

0.  979 

0.  999 

1.019 

1.035 

1.046 

1.062 

0.  95 

17.6 

15.  9 

14.  4 

13.  1 

12.  2 

11.  7 

10.  9 

8.  5 

Units:  CK),  M^/Q  Ckg/Kw),  L  Cm).  D  Cm),  a  Cmicrons) 


LlI 


Results  for  100 
CBeta  = 


Realistic  Values 

=  100.0) 


Mission 

Time 


10  years 


20  years 


30  years 


Units: 

*  poss: 


p 

ump  Specific 

:  Mass 

Ckg/ Ckg/ sec) ) 

tan 

15.0 

20 . 0 

25.0 

30.0 

35.0 

40.0 

250 

246 

243 

240 

237 

235 

232 

0.  828 

0.  886 

0.  939 

0.  990 

1.039 

1.086 

1.  131 

218 

222 

225 

228 

231 

233 

237 

0.  87 

0.89 

0.  90 

0.  91 

0.  93 

0.  93 

0.  95 

32.7 

28.  6 

26.0 

23.  8 

21.  8 

20.  6 

19.0 

249 

245 

242 

239 

236 

234 

232 

0.  846 

0.  902 

0. 955 

1.005 

1.053 

1.099 

1.  144 

219 

223 

226 

229 

232 

235 

237 

0.  88 

0.  89 

0.  90 

0.  92 

0.  93 

0.  94 

0.  95 

31.  6 

27.  7 

25.  2 

23.  1 

21.  2 

20.  1 

19.0 

247 

244 

241 

238 

236 

233 

231 

0.  862 

0.  917 

0.  969 

1.019 

| 

1.066 

1.  112 

1.  156 

221 

224 

227 

230 

232 

236 

238 

0.  89 

0.  90 

0.  91 

0.  92 

0.  93 

0.  94 

0.  95 

29.  5 

26.  8 

24.  5 

22.  4 

21.  2 

19.  5 

18.  5 

246 

243 

240 

237 

235 

233 

231 

0.  879 

C.  933 

0.  984 

1.033 

1.080 

1.  125 

1.  169 

222 

225 

228 

231 

233 

236 

238 

0.  89 

0. 90 

0.  91 

0.  93 

0.  93 

0.  94 

0.  95 

28.  6 

26.0 

23.  8 

21.  8 

20.  6 

19.  5 

18.  5 

T?  Cl 
Lble  t: 

O,  Hp, 

vpo  in 

'fit  Ckg 
Knaoo 

/Kw),  1 
study 

-  Cm), 

D  Cm), 

a  Cm 

Knapp 


246 
0.  979 
221 
0.  88 
12.0 


2 


0.  995 


224 
0.  90 
11.  4 


24: 

1.01' 

22 

*0.09 

11. 


TABLE  VI 


Results  for  1000  Kw  Using  Optimistic  Values 
(Beta  =  0.  10,  me  =  40.0) 


Mission 

Pump  Specific  Mass 

Ckg/ (kg/sec) ) 

Time 

mm 

mm 

20.  0 

ggHjfl 

mm 

40.  C 

Knapp 

o  t2= 

219. 

216 

213 

211 

208 

206 

204 

222 

Mj/Qs 

0.  570 

0.  606 

0.  641 

0.  675 

0.  708 

0.  740 

0.  772 

0.  528 

L= 

797 

810 

822 

831 

844 

854 

863 

781 

D= 

3.  19 

3.  24 

3.  29 

3.  32 

3.  38 

3.  41 

3.  45 

3.  12 

a= 

43.  4 

40.  5 

37.  8 

36.  2 

33.  9 

32.  4 

31.  1 

24.0 

10  years 

218 

215 

212 

210 

208 

2C6 

2C4 

222 

0.  580 

0.  616 

0.  651 

0.  685 

0.  717 

0.  749 

0.  781 

0.  539 

801 

814 

827 

835 

844 

854 

863 

781 

3.21 

3.  26 

3.  31 

3.  34 

3.  38 

3.  41 

3.  45 

3.  12 

42.  4 

39.  6 

37.0 

35.  4 

33.  9 

32.  4 

31.  1 

24.0 

20  years 

217 

214 

212 

209 

207 

205 

203 

222 

0.  591 

0.  626 

0.  661 

0.  694 

0.  727 

0.  759 

0.  790 

0.  550 

805 

818 

827 

840 

849 

858 

068 

781 

3.  22 

3.  27 

3.  31 

3.  36 

3.  40 

3.  43 

3.  47 

3.  12 

41.  4 

38.  7 

37.  0 

34.  6 

33.  1 

31.  7 

30.  4 

24.0 

30  years 

216 

213 

211 

209 

206 

2C4 

203 

219 

0.  601 

0.  637 

0.  671 

0.  704 

0.  736 

0.  768 

0.  799 

0.  561 

810 

822 

831 

840 

854 

863 

868 

793 

3.  24 

3.  29 

3.  32 

3.  36 

3.  41 

3.  45 

3.  47 

3.  17 

40.  5 

37.  8 

36.  2 

34.  6 

32.  4 

31.  1 

30.  4 

22.  7 

Units:  T?  CK),  Mq,/Q  Ckg/Kw),  L  Cm),  "  Cm),  a  (microns) 
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TABLE  VII 


Results  for  1000  Kw  Using  Realistic  Values 
(Beta  =  0.20,  me  -  100.0) 


Mission 

Pump  Specific  Mass 

Ckg/Ckg/sec) ) 

Time 

man 

wsm 

nnu 

B» 

Knapp 

o  t2= 

236 

234 

232 

230 

228 

226 

225 

237 

mt/Ui= 

1.025 

1.071 

1.  115 

1.  158 

1.  200 

1.  242 

1.  282 

1.032 

L= 

735 

742 

749 

756 

763 

770 

774 

727 

D= 

2.  94 

2.  97 

2.  99 

3.02 

3.05 

3.08 

3.09 

2.  91 

a= 

67.  1 

63.  4 

60.  1 

57.0 

54.  1 

51.5 

50.  2 

32.0 

10  years 

235 

233 

231 

229 

228 

226 

224 

237 

1.038 

1.084 

1.  128 

1.  171 

1.  213 

1.  253 

1.  293 

1.  046 

738 

745 

752 

759 

763 

770 

778 

727 

2.  95 

2.  98 

3.01 

3.04 

3.05 

3.08 

3.  11 

2.  91 

_ _ 

65.2 

61.7 

58.  5 

05.5 

54.  1 

51.5 

49.0 

32.0 

20  years 

235 

233 

231 

229 

227 

225 

224 

237 

1.052 

1.097 

1.  140 

1.  183 

1.  225 

1.  265 

1.305 

1.068 

738 

745 

752 

759 

767 

774 

779 

727 

2.  95 

2.  98 

3.00 

3.04 

3.07 

3.  10 

3.  11 

2.  91 

65.  2 

61.  7 

58.  5 

55.  5 

52.  8 

50.  2 

48.  9 

32.0 

30  years 

234 

232 

230 

228 

227 

225 

223 

234 

1.065 

1.  110 

1.  153 

1.  195 

1.  237 

1.  277 

1.  316 

1.074 

742 

749 

756 

763 

767 

774 

782 

738 

2.  97 

2.  99 

3.02 

3.05 

3.07 

3.  10 

3.  13 

2.  95 

63.  4 

60.  1 

57.0 

54.  1 

52.  7 

50.  2 

47.  8 

30.0 

Units:  CK),  ty^/Q  Ckg/Kw),  L  Cm),  D  Cm),  a  Cmlcrons) 
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TABLE  VIII 


Results  for  10000  Kw  Using  Optimistic'  Values 
(Beta  =  0.10,  me  =  40.0) 


Mission 


Pump  Specific  Mass  (kg/ (kg/sec)) 


Time 

KUUKUUBSBil 

30.0 

35.0 

40.0 

Knapp 

o 

H 

ro 

II 

198 

197 

196 

194 

193 

192 

191 

201 

Mp/Q= 

0.  901 

0.  930 

0.  959 

0.  987 

1.016 

1.043 

1.071 

0.  845 

L= 

2823 

2839 

2855 

2889 

2906 

2923 

2940 

2759 

D= 

11/3 

11.  4 

11.  4 

11.  5 

11.  6 

11.  7 

11.  8 

11.0 

a= 

86.  8 

85.0 

83.  3 

80.0 

78.  4 

76.  9 

75.  4 

53.0 

10  years 

198 

197 

195 

194 

193 

192 

191 

201 

0.  909 

0.  938 

0.  967 

0.  996 

1.024 

1.052 

1.079 

0.  857 

2823 

2839 

2872 

2889 

2906 

2923 

2941 

2759 

11.  3 

11.  4 

11.  5 

11.  6 

11.  6 

11.  7 

11.  8 

11.0 

86.  8 

85.0 

81.  6 

80.0 

78.  4 

76.  9 

75.  4 

53.0 

20  years 

197 

196 

195 

194 

192 

191 

190 

201 

0,  918 

0.  947 

0.  976 

1.004 

1.032 

1.059 

1.087 

0.  866 

2839 

2855 

2872 

2889 

2923 

2941 

2959 

2759 

11.  4 

11.  4 

11.  5 

11.  6 

11.  7 

11.  8 

11.  8 

11.0 

85.0 

83.  3 

81.  6 

80.  0 

76.  9 

75.4 

73.  9 

53.0 

30  years 

197 

196 

194 

193 

192 

191 

190 

198 

0.  926 

0.955 

0.  984 

1.012 

1.040 

1.068 

1.095 

0.  875 

2839 

2855 

2889 

2906 

2923 

2941 

2959 

2807 

11.  4 

11.  4 

11.  6 

11.  6 

11.  7 

11.  8 

11.  8 

11.  2 

i 

85.0 

83.  3 

80.  0 

78.  4 

76.  9 

75.  4 

73.  9 

51.0 

Units:  T2  CK) ,  M^/Q  Ckg/Kw),  L  Cm),  D  Cm),  a  Cmicrons) 
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TABLE  IX 


Results  for  10000  Kw  Using  Realistic .Values 
(Seta  =  0.20,  me  =  100.0) 


.Mission 

Pump  Specific  Mass 

C kg/ C kg/ sec) ) 

■■■ 

Time 

IEBS1 

iHia 

j&am 

i£B*l 

KCTP1 

CESS! 

40.0 

BBSS 

it 

C\J 

H 

o 

216 

215 

214 

213 

212 

211 

210 

*232 

MT/<i= 

1.48 

1.52 

1.55 

1.  59 

1.  62 

1.  65 

1.69 

1.  44 

L= 

2560 

2573 

2587 

2600 

2614 

2628 

2642 

2507 

D= 

10.  2 

10.  3 

10.3 

10.  4 

10.5 

10.5 

10.  6 

10.0 

a- 

128 

125 

122 

119 

117 

114 

112 

72.0 

10  years 

216 

215 

214 

213 

212 

211 

210 

216 

1.49 

1.53 

1.56 

1.59 

1.  63 

1.  66 

1.  70 

1.  45 

2560 

2573 

2587 

2600 

2614 

2628 

2642 

2546 

10.  2 

10.  3 

10.3 

10.  4 

10.  5 

10.  5 

10.  6 

10.  2 

128 

125 

122 

119 

117 

114 

112 

68.0 

20  years 

216 

215 

214 

213 

212 

211 

210 

216 

1.50 

1.  54 

1.57 

1.61 

1.  64 

1.  67 

1.  71 

1.  46 

2560 

2573 

2587 

2600 

2614 

2628 

2642 

2546 

• 

10.  2 

10.  3 

10.3 

10.  4 

10.5 

10.  5 

10.  6 

10.  2 

128 

125 

122 

119 

117 

114 

112 

68.0 

30  years 

215 

214 

213 

212 

212 

211 

210 

216 

1.51 

1.  55 

1.58 

1.  62 

1.  65 

1.  68 

1.  72 

1.  47 

2573 

2587 

2600 

2614 

2614 

2628 

2642 

2546 

10.  3 

10.  3 

10.  4 

10.  5 

10.  5 

10.  5 

10.  6 

10.  2 

125 

122 

119 

117 

117 

114 

112 

68.0 

Units:  T„  CK),  typ/Q  Ckg/Kw),  L  Cm),  D  Cm),  a  Cmicrons) 
*  possible  typo  in  Knapp  study 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  study  was  a  parametric  analysis  of  the  Liquid 
Droplet  Radiator  CLDR)  using  a  new  pump  specific  mass  (mp) 
term  defined  as  pump  mass  per  liquid  mass  flow  rate  to 
calculate  the  minimum  LDR  system  mass  per  heat  radiated  as  a 
function  of  the  average  droplet  temperature  at  the 
collector.  The  characteristics  of  total  system  specific 
mass,  stream  length,  stream  diameter  and  droplet  radius  were 
calculated  as  a  function  of  the  average  droplet  temperature 
at  the  collector  for  the  minimum  system  mass  for  a  range  of 
mp  values  from  10.0  to  40.0  kg/Ckg/sec)  by  increments  of 
five.  The  results  of  this  study  indicate  the  new  specific 
pump  mass  term  provides  a  physically  meaningful  term  for 
design  engineers  to  use  on  future  LDR  systems  by  offering  a 
range  of  pump  masses  depending  on  pressure  loss  and  flow 
rate  design  requirements. 

The  assumed  values  of  the  parameters  used  in  this  study 
should  be  the  subject  of  futher  research  to  provide  possible 
improvements  to  the  values.  These  parameter  values  include 
the  ratio  of  liquid  mass  in  the  reheating  station  to  droplet 
mass  in  the  stream,  the  average  emi3sivity  of  the  droplets, 
the  ratio  of  droplet  kinetic  energy  to  heat  rejected,  the 
conservative  black  body  view  factor  for  a  droplet  at  stream 
center,  the  change  of  vapor  pressure  of  the  droplets  with 
changing  temperature,  and  the  specific  mass  term  for  the 
generator  and  collector.  A  previous  study  by  Knapp  Cl) 
using  the  3ame  values  for  these  parameters  with  a  different 
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pump  specific  mass  defined  as  kg-sec  /m  provided  results 
within  a  few  percent  agreement  with  the  characteristic 
results  of  this  study  for  total  system  specific  mass, 
average  droplet  temperature  at  the  collector,'  and  LDR  stream 
length  and  diameter  calculated  with  a  mp  value  of  10.0 
kg/Ckg/secD  primarily  for  the  heat  loss  rates  of  100,  1000, 
and  10000  Kilowatts  in  both  the  optimistic  and  realistic 
cases.  Knapp's  results  for  these  characteristics  at  10 
Kilowatts  best  fit  the  values  obtained  with  mp  =  10.0  for 
the  optimistic  case  and  best  fit  the  characteristic  values 
obtained  with  mp  =  40.0  for  the  realistic  case.  Also, 
Knapp's  droplet  radius  was  low  in  all  cases  by  approximately 
one  half  the  calculated  values  of  this  study. 

Additionally,  the  total  system  specific  masses 
calculated  in  this  study  were  lower  than  the  10  kg/Kw  value 
C 1 : 14)  for  heat  pipe  or  fluid  loop  radiators  at  the  300  K 
temperature  level.  This  reduced  mass  per  heat  radiated 
makes  the  LDR  superior  to  traditional  heat  pipe  and  fluid 
loop  radiators  for  high  power  space  applications. 

Further  analysis  is  needed  to  understand  differences  in 
the  results  at  the  10  Kilowatt  heat  loss  rate.  Additional 
analysis  at  different  heat  loss  rates  may  provide  some 
insight.  Also  higher  droplet  generator  temperatures  using  a 
liquid  metal  might  be  used  to  observe  the  effect  on  the 
calculated  characteristics.  Continued  work  in  each  of  these 
areas  will  provide  the  necessary  research  to  enable  the  LDR 
to  become  the  radiator  of  choice  for  future  space  missions. 
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Appendix  A 


Computer  Program  Glossary  and  Listing 


50  * ••••••••ssssesseoLOSSAHTee ••*••••■••••••••••••• •••••••• 

60  • 

70  '  A  -  radius  of  spherical  particles 

60  '  AH  -  variable  in  quadratic  foraula 

90  '  B  -  variable  (function  of  ETA  and  EPS) 

100  '  BE  -  variable  in  quadratic  foraula 

110  '  BETA  -  rttlo,  liquid  aaaa  in  reheating  atation  to  mass  in  atreaa 

120  '  C  -  apeoifio  heat  of  dropleta 
130  *  CE  -  variable  in  quadratic  foraula 

1 A0  '  0  -  diaaeter  of  droplet  atreaa 

150  •  EPS  -  eaisslvlty  of  dropleta 

160  '  ETA  -  black-body  view  faotor  for  a  droplet  at  atreaa  center 

170  *  P  -  gray-body  view  factor  for  a  droplet  at  atreaa  oenter 

180  '  QAHA  -  ratio  of  klnetio  energy  to  the  radiated  energy 

190  '  L  -  length  of  droplet  atreaa'  along  s-axls 

200  •  HPQE  -  aaaa  per  watt  of  radiated  heat  losa 

210  •  PI  -  J.1A159 

220  •  PT  *  vapor  praaaure  of  droplet  at  teaperature  T1D 

230  *  QT3D  -  total  rate  of  heat  loaa  to  apace  froa  droplet  atreaa 

2*0  •  P  -  gaa  constant 

250  ’  EH0  -  density  of  droplet 

260  ’  S10  -  Stefan-Boltzaann  constant,  5.67E-GB  vatts/aeter*2-kelvin*A 

270  *  SH0C  -  specific  asss  of  droplet  generator  and  oollector 

280  •  3MP  -  specific  aass  of  puaps 

290  *  TA0  -  alssion  lifatiae 

300  *  T1D  -  initial  teaperature  of  droplets 

310  ’  T2D  -  final  teaperature  of  dropleta  at  end  of  streaa 

320  *  TV  -  teaperature  used  in  fitting  vapor  pressure  curve 

330  •  V  -  velocity  of  the  streaa 

3*0  ;  W1-W5  -  variables  in  HPQE  equation 

350  * 


10  •  PROGRAM;  TSMYLDR  -OSES  A  NEW  TERM  FOR  TKE  PUMP  MASS  IN  MPQR 
15  •  EQUATION.  FART  LIRE  580  FOR  T2D  TO  GET  MINIMUM  MPQR. 

20  >  QTSD,  TAU,  AND  SMP  MUST  ALSO  BE  DESIGNATED  AT  INPUT. 

30  •  THE  VALUES  FOR  BETA  AND  SMGC  MOST  BE  CHECKED  (LINES  RAO  A  *60) 
35  '  FOR  EITHER  THE ‘OPTIMISTIC  VALUE'  OR  THE  REALISTIC  VALUES. 

AO  • 

350  •••■••DEFINE  VARIABLES  FOR  DON  705  SILICON  OIL  FLUID 


•KELVIN 
•KELVIN 

•JOULE/KILOQRAM-KBLVIN 
•KILOGRAM/CUBIC  METER 
•ASSUMPTION 
•ASSUME 
•ESTIMATED 

•ASSUME  (KILOGRAMS/SQUARE  METER) 

•  JOULE/KELVIN 
•NEWTONS/SQUARE  METER 

•  WATTS /ME7 ER*2-KELVIH“A 

510  • 

520  •••#Mti*»»«»i««i»i«»t»iV»»i»««#PMFORM  CALCULATIONS 
530  • 

5A0  CLS 

550  GOSUB  820  'CALCULATE  2TA  AND  F 

560  OOSUB  690  'ASK  FOR  INPUT 

570  I*  1+1  :  J*  J+1  'INITIALIZE  COUNTERS 

580  T2D«  230  +  I  'VAR I  T2D  TO  FIND  MINIMUM  MPQR 

590  PRINT $PRINT"T2D*;T2D 

600  E*  (8*(1-ETA))/((1/ETA)+(1/EPS)-1)  :'PRINT"B»B 
610  GOSUB  950  'CALCULATE  MPQR 

620  IF  LD  >  250  THEN  LD.250  :  OOSUB  1180  'CALC  MPQR  FOR  L/D*  250 

630  IF  MPQR(J-I)*  0  THEN  MPQR(J-I)*  9.999999E-21 

6A0  IF  ABS((HPQA(J)-MPQR(J-1))/MFQRU-1))  >  .0001  THEN  GOTO  570  'TO  CONVERGE 
650  QTSD«(((8»(1-ETA)»D*L»SI0»T1D*A)/((1/ETA)+(1/EPS)-1))»((3»(1-(T2D/T1D)))/ 

( (T2D/T1D)*-3-1 > ) )  sPRIMT"L650  QTSD" ;QTSD 


660  LPRINT"  CONVERGED  TO  THESE  VALUES" 

670  END 
680  • 

690  « *••»•••»»»•••»•*•••••••«•* request  INPUT  PARAMETERS 

700  • 


710  INPUT*QTOT  IN  VATTS" ;QTSD : PRINT 
720  INPUT"TIME  IN  TEARS" | TAU : PRINT 
730  INPUT»PUMP  SPECIFIC  MASS  (KG  PER  KO/SEC)"jSMP 
7*0  LPRINT"  Q  TOT*" ; QTSD ; "HATTS" SPRINT 

750  PRINT"Q  TOT*" ; QTSD; "WATTS" SPRINT 

760  LPRINT"  TIME* ";TAU; "YEARS  SMP*» JSMP ; "KG  PER  KG/SEC" 

770  PRINT"TIME»»jTAU; "TEARS* SPRINT 

780  TAU*  TAU«3. 1536E+07  'CONVERT  TIME  IN  TEARS  TO  SECONDS 

790  LPRINT*  T2  MT/Q  L  D  A" 

800  RETURN 
810  ' 


380  PI*  MATN(1) 
390  T1D*  300 
*00  TV*  15000 
RIO  C*  1670 
*20  RHO*  1000 
*30  GAMA*  .005 
**0  BETA*  .1 
*50  EPS*  .7 
*60  SMGC*  *01 
*70  R*  15.22 
*80  PT*  3.001E-08 
*90  SIG*  5 . 67E-08 
500  DIM  MPQR(IOO) 
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CALC  ETA  AND  P 


8?0  'HHHIHIHHIHHHeilHIHHI 

830  • 

SAC  IP  EPSn  0  THEN  EfSc  1.00001E-30 
850  IP  EPS*  1  THEN  ETA*  .5  :  GOTO  900 
860  AH*  ( ( 1 l/EPS)-1 I ) 

870  BE*  21 
880  CB.  >11 

890  ETA*  (-BE  ♦  SQR(BE*2  -  A*AH«CE) )/(2«AH)  I 

900  P«  1/<<1/*PS)+<1/STA>-1)  ! 

910  PRINT  USINC’EKS *##.##  F> ##.###«  ETA. M .«»*#*; EPS ;F;ETA 

920  PRINT 

930  RETURN 

5A0  • 

950  •••••••••••••••••CALC  HASS  PER  NATT  HEAT  LOSS-  MPQR 

96C  *  ; 

970  V*  SQR(0AMA*2*C*(T1D-T2D) )  :  PRINT" ¥•» ; "0AMMA"0AHA  J 

980  VI.  (3*QISD'(1/3))/(2»(1-(T2D/T1D))*(A/3))  s'PRINT"V1"W1 

990  V2*  (<(UBETA)»<((T2D/T1D)*-3>-1>>/<B»C*¥»T1D»3IQ»T1D*«l)>*(2/3>:»PRINT  V2  ! 
1000  V3*  ( (SHOC»PI)/ 18 )*( 1/3)  J'PRINT*V3*H3  i 

1010  V*«  3MP/(C»(T1D-T2D))  ! ' PRINT*V4*VA  ! 

1020  V5*  (ETA»PT»(T1D/T¥)«TAU)/(F«3IQ»T1D**»SQR(2»PI»R»T1D)»(1-(T2D/T1D)>) 

1030  • PRINT"V5"W5 

10M0  MPQR(K).  (H1*H2*H3)«V«-»V5 

1050  PRINT  USING*MFQR*#f .## f ##  KILOORAMS/KILOHATT  EQN  23" ;KPQR(K)» 1000 
1060  L.  ((PI»C«V»T1D»QTSD»((T2D/T1D)*-3-1)*2,SM0C)/(18»B*2»(UBETA)» 
(SIO«T1D*«)*2»(1-(T2D/T1D))))*(1/3) 

1070  D.  QTSD“(2/3)»(((2»(UBETA))/(3#B*C»¥»T1D»SHOC»SIO»T1D‘A))“(1/3))» 
((((T2D/T1D)*(l/3))/( 1-(T2D/T1D) )“ ( 2/3) ) ) 

1080  A.(9»P/RHO)»(1/(C»t»T1D»B)“(2/’3))#((M,SHOC»SIO»T1D*A)/(18»(nBE-A)))',(1/3) 
•QTSD*(1/3)/(( (T2D/T 1D)"-3-1 )*( 1/3)*(( 1-(T2D/T1D) )“( 1/3))) 

1390  LD*  L/D 

1100  'PRINT  *L*" jL; "HETERS" 

1110  'PRINT’D." ;D {"METERS" 

1120  'PRINT"A.";A;BHETERS" 

1130  ' LPRINT  SPC(105  T2D;  HPQR(R)» 1000;  L;  D;  A 
1*.A0  PRINT  T2D;  MPQROO*  1000 ;  L;  D;  A 
1150  PRINT’EQN  13  L/D«";LD  SPRINT 
1160  RETURN 
1170  ' 


1180  •••■•••••••••••CALC  MASS  PER  WATT  BEAT  LOSS-  MPQR  (POR  L/P«  250) 

1190  ' 

1200  T«  SQR(0AMA*2»C»(T10-T2D) )  s  PRIRT "7"T; "GANN A "QAM A 
1210  Vie  (<WRETA)*3QRUB)»<(T2B/T1»D*-3-1)*(1/2)«QTSB*<1/2))/ 
(C*I*T1P#SQB(3*R,3I0*T1B“A)*{ 1-(T2B/T1B) )*( 3/2)) 

1220  »PRIRT"1M"I  I 

1230  M2»  (PI*SMOC*( (T2B/T1B)“-3-1 ) )/( 12»§»SI0*T1D*R»(LD),( 1-(T2B/T1P) ) ) 
1280  • PRIBT*N2*V2 

1250  M3«  SMP/(C»<T10-T2D))  : •PRIHT"B3"B3 

1260  «*«  \S?A,Pr»(T10/TT)*TA0)/{P»SI0*T10*R»SQR{2»PI»R,T10)#( 1-CT2D/T1D) ) ) 

1270  •PRIBT"H««VA 

1280  MPQR(J)«  HWM2«N3*N* 

1290  PRIRT  03IHG"MPQRs  ##.####  KILOGRAMS/ KILOWATT  EQR  26";M?a5{ if  1000 

1300  D.  (QIS0'‘(1/2)«(<T2D/T1B)*-3-1)'(1/2))/((3#»,LB*3I5SirlB*R)*{1/2)* 

( 1-(T2B/T1B))“( 1/2) ) 

1310  L*  B«LB 

1320  An  (9,SIO#P#L/(C*RHO*y)),{T1B*3/( (T2B/T1B)*-3-1 )) 

1330  ' PRIBT"L«*iL;aHKTERS" 

1380  •PSIBT*P«";B; "METERS" 

1350  • PRIBT"A«*;A;"METKRS" 

1360  LPRIRT  SPC(IO)  T2D;  MPQR( J)» 1000 ;  L;  0;  A 
1370  PRIRT  T2B;  MPQR(J)*1000;  L;  B;  A 
1380  PRIRT-EQR  18  L/B";L/B 
1390  RETURH 
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■  Appendix  B 


Pump  Specif ic  Mass 

The  work  in  Reference  1  used  two  values  for  pump 

2/3  2 

specific  mass  with  units  of  kg-sec  /m  .  The  optimistic 

2/7  2 

value  was  given  as  1000  kg-sec  '  J/m  and  the  realistic  value 
2/3  2 

was  2000  kg-sec  "7m  .  These  values  made  little  physical 
sense  from  a  design  point  of  view.  This  study  used  a  new 
pump  specific  mass  (mp)  for  a  centrifugal  pump  in  terms  of 
the  pump  mass  per  mass  flow  rate.  The  new  mp  values  were 
varied  from  10.0  to  40.0  as  shown  in  Table  I.  This  range  of 
values  appeared  reasonable  from  available  stainless  steel 
pump  data  (8:583’).  Additionally,  Mattick  and  Taussig  (4:14) 
describe  a  droplet  radiator  system  where  the  value  of  mp  was 
found  to  be  10.05  kg/(kg/sec).  Therefore  a  range  of  mp 
values  was  used  in  this  study  to  observe  the  effect  on  the 
operating  characteristics  of  the  LDR.  A  sample  of  mp  values 
shown  below  was  calculated  from  Reference  8. 


TABLE  X 

Pump  Specific  Mass  Sample  Values 


Pump  Model 

mp  (kg./ (kg/ sec) ) 

J7005-30 

9.  12 

J7005-35 

10.07 

J7C05-20 

13.  26 

J7005-10 

18.  87 
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available  industry  standard  pumps.  The  average  droplet 
temperature  at  the  collector  was  calculated  to  minimize  the 
LDR  system  mass  for  heat  loss  rates  of  10  Kw,  100  Kw, 

1000  Kw,  and  10,000  Kw  for  mission  lifetimes  of  zero,  ten, 
twenty,  and  thirty  years.  The  initial  droplet  temperature 
at  the  generator  was  fixed  at  300  degrees  Kelvin.  A  silicon 
oil,  Trimethylpentaphenyltrisiloxane  (DOW  705),  was  modeled 
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